A New Quadruply Lensed Quasar: SDSS J125107.57+293540.5 
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We report the discovery of a quadruply imaged quasar, 
SDSS J125107.57+293540.5, selected from the Sloan Digital Sky Survey. 
Follow-up imaging reveals that the system consists of four blue point-like 
components in a typical cusp lens geometry surrounding a central red galaxy. 
^ . The source redshift is 0.802 and the lens redshift is 0.410. The maximum 

image separation between the lensed components is 1'.'79. While the image 
configuration is well reproduced by standard mass models with reasonable 
. parameter values, the flux ratios predicted by these models differ from the 
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observed ratios in all bands. This is suggestive of small-scale structures in this 
lens, although the definitive identification of the anomaly requires more accurate 
photometry and astrometry. 

Subject headings: gravitational lensing — quasars: individual (SDSS J125107.57+293540.5) 



Introduction 



Nearl y one hundred gravitationally lensed quasars have been found to date (see the 
review by iKochanekl l2006l ) , and they have become a unique astronomical and cosmological 
tool. Their usefulness is particularly enhanced by constructing a statistical sample with a 
well- understood selection fun ction with the hel p of extensive and homogeneous surveys such 
as the HST Snapshot Survev (IMaoz et al.lll993l) or the Cosmic Lens All Sky Survey (CLASS; 
Myers et al.ll2003l : iBrowne et al.ll2003l ). However, the number of strongly lensed quasars in 
systematic surveys is still limited. To construct a l arger lens sample, we are conducting 
the SPS S Quasar Lens Search (SQL S; lOguri et al.l l2006l ) based on spectroscopic quasar 
catalogs (ISchneider et al.l l2005l . 120071 ) produced by the Sloan Digital Sky Survey (SDSS; 
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20011 ). A significant fraction of our lens candidates still require follow-up observations, so we 
expect to discover more lensed quasars in the coming years. 

In this paper, we report the discovery of a quadruply lensed quasar, SDSS J125107. 57+293540. 5 
(SDSS J1251+2935), in the course of the SQLS. Quadruple lenses, which constitute roughly 
one-third of all lensed quasars, are not only visually interesting objects due to their charac- 
teristic morphologies but also scientifically important and useful objects. Four image lenses 
provide many more constraints on the mass distribution of the lens. In particular, the im- 
age fiux ratios expected for smooth, central potentials are relatively well-defined and the 
differences between the observed and model flux ratios can be used to study gravitational 



substructures (e.g., iMao fc Schneider I Il998l ). Anomalous flux ratios, where the observed 



and model image fluxes are in signiflcant disagreement are fairly common and are due to the 
combined effects of the stars in the lens galaxy and satellite halos. In radio lenses this has 
been used to argue that the abundance of satellite halos is consistent with the predictions o f 



standard CDM models fe.g.. iMetcalf fc Madaull200ll : IChiball2002l : iKochanek fc Dalalll2004[ ) 



while in optical lenses the effects of substructure and microlensing by the stars must be 
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disentangled (e.g. 



Chiba et al. 


2005; 


Keeton et al. 


2006; 


Morgan et al. 


2006) 



In the next section, we briefly describe our algorithm for candidate selection from the 
SDSS data. Details of follow-up imaging using the University of Hawaii 2.2-meter (UH88) 
telescope are provided in § [31 We investigate mass models of the system in § HI and 
summarize our results in § [5l Throughout the paper we assume a cosmological model 
with matter density Qm = 0.27, cosmological consta nt = 0.73, and Hubble constant 
h = i7o/(100km s-^Mpc"^) = 0.7 JSpergel et al.ll2003h . 



2. SDSS Data 



The SDSS is conducting both a photorn etric survey (iGunn et al.l Il998l : iLupton et al. 

19991 : Istoughton eraDl2002l : iTucker et al]l2006h in flve broad-band optical fllters (IFukugita et al, 
19961 ) and a spectroscopic survey with a multi-flber spectrograph covering 3800 A to 9200 A 
at a resolution of R ~ 1850 . The SDSS uses a dedicated wide-fleld (3° fleld of view) 2.5-m 



telescope (iGunn et al.ll2006l ) at the Apache Point Observatory in New Mexico, USA, covering 
10,000 square degrees of the sky approximately centered on the North Gala ctic Cap. The 
imaging data are processed by the photometric pipeline (ILupton et al.ll200ll ). and spectro- 
sc opic quasar targ e ts are selected from the imaging data according to the algorithm described 
by [Richards et al. J2002[l_ Fibers for t he sp ectroscopic observations are assigned according to 
the tiling algorithm of Isianton et al. ( 2003). The imaging data have an astrometric accuracy 
better than about O'.'l rms per coordinate (jPier et al.ll2003l ) an d photometric zeropoint errors 



less than about Ci 



2002 



data; 



Ivezic et al. 



.03 rn agnitude over the entire survey area (IHogg et al.l l200ll : ISmith et al. 



2004r). SDSS J1251-I-2935 is contained in Data Release 5 (for the imaging 



Adelman-McCarthyl 120071 ) and later (for the spectroscopic data). 



SDSS J125H-2935 was sele cted as a lensed q uasar candidate using the morphological 
selection algorithm described in lOguri et al.l ( 20061 ). The algorithm uses the SDSS morpho- 
logical classiflcation parameter objc_type and likelihood star_L that an object is fltted by a 
point-spread function (PSF). Although lensed quasar systems with small image separations 
are classifled as single objects in the SDSS data, the proflles are extended and are not con- 
sistent with either PSF proflles or simple galaxy proflles. Therefore, small separation lensed 
quasar candidates are selected as objects that have very small values of star_L. In addition, 
SDSS J12 51-I-2935 satisfles additional selection requirements based on flts to the image with 
GALFIT ( jPeng et al.l 120021 ) which are applied to exclude false positives by single quasars. 
While some single quasars can pass the initial selection step, flts to such systems using two 
PSFs lead to either very large magnitude differences or very small image sepa rations that 
are indicative of systematic errors rather than a gravitational lens (see § 5 of lOguri et al. 
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(120061)). 

The SDSS z-band image of SDSS J1251+2935 is shown in Figuredl The PSF magnitudes 
of SDSS J1251+2935 (after correcting for Galactic extinction) are 19.83 ±0.04, 19.38 ±0.02, 
19.13±0.03, 18.85±0.03, and 18.42±0.04 in u, g, r, i, and z, respectively. SDSS J1251±2935 
is spectroscopically confirmed as a quasar at z = 0.802 (see Figure[2]for the SDSS spectrum). 
In this spectrum, we can also see a series of absorption lines that indicate the presence of a 
bright early-type galaxy at z = 0.410 ± 0.001. The presence of both components is a strong 
indication that this system is a gravitational lens. 



3. Imaging Follow-up Observations 

Deeper and higher- resolution optical images of SDSS J1251±2935 were obtained on 2006 
April 25 (0'.'8 seeing) and 2006 May 3 (I'.'O seeing) using the 8k mosaic CCD camera (UH8k, 
pixel scale of 0'.'232 pixel"^) and the Orthogonal Parallel Transfer Imaging Camera (Optic, 
pixel scale of 0'.'137 pixel"^) at the UH88 telescope, respectively. We took V- and /-band 
images (270 sec exposure for each band) with the UH8k, and B-, R-, and J-band images 
(400 sec exposure for each band) with the Optic. Because the night of UH8k imaging was 
not photometric, we did not photometrically calibrate the UH8k images. We binned (2 x 2) 
the Optic images and used them for the astrometry and photometry of the system. The 
2x2 binned Optic images (BRI) and the original UH8k image {V) are shown in the left 
panel of Figure [3] and the upper panels of Figure HI 

We used GALFIT to model these images with a series of models of increasing complexity. 
The only model that works well consists of 4 point sources and a central galaxy modeled 
with a Sersic profile. In Figure [3], we demonstrate that subtracting 4 fitted PSFs from the 
Optic /-band image leaves an extended object or vice versa. If we further subtract the galaxy 
component, there remains virtually no residual. The lower panels of Figure H] summarize the 
residuals after subtracting the best models for each band's image. The galaxy flux is well- 
determined only in the R and /-band images, so we neglected this component in the B and 
V^-band fits even though there are hints of its presence in the l^-band residuals. We label the 
4 point sources A-D in order of increasing /-band magnitudes and the central galaxy as G. 
We estimated that the galaxy has an effective radius of I'.'ll ± 0'.'38, ellipticity of 0.28 ±0.09, 
the Sersic index of 2.4±1.1, and a major axis position angle of 26°±5°. The R — I color of G, 



0.89 ± 0.30, is consistent with an early-type galaxy at the measured redshift (iFukugita et al. 



19951 ). Table [T] summarizes the relative astrometry and photometry of the system, where we 
defined the errors from the scatter between the fits using 6 different PSF templates rather 
than the smaller statistical uncertainties of the individual fits. The errors are an order of 
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magnitude larger than other systems where images are resolved even by visual inspection. 
We plot the color-color diagram {B — R and R — I) of the 4 point-like components in Figure [5] 
and the flux ratios between the point- like components in Figure O Although the colors (and 
flux ratios) of the lensed images have large scatter among the images, they are consistent 
with each other given the large uncertainties. 

In Figure [H we also labeled the nea rby galaxies by their R — I colors as measured 
using SExtractor (IBertin fc Arnoutsi Il996l ). There are several galaxies with colors similar 
to the lens galaxy to the Southwest (1.0 < R — I < 1.1), and a larger group of generally 
bluer galaxies (0.6 < R — I < 1.0) ~ 40" to the North. This sugge sts that the lens is 
associated with a group , as is quite common among lensed quasars (e.g.. iFassnacht &: Lubin 
2OO2I : lOguri et al.lboosl : lOgurilbood : Iwilhams et aUbooeh . 



4. Mass Modeling 



We modeled the system using the /ensmorfe/ package (lKeetonll200ll ) to determine whether 
reasonable mass distributions can reproduce the observations. We first used the 7 parameter 
singular isothermal ellipsoid model (SIE; the Einstein radius Re, the ellipticity e and its 
position angle 9e, and the positions of the galaxy and the source quasar) to fit the relative 
positions of the A-D and G components measured from the i?-band image. The model fits 
the data well with Xrcd = X^/ dof = 1.1 for dof = 3 degrees of freedom. The model parameters 
are presented in Table [2] and the fit is illustrated in Figure [71 The predicted position angle 
of the lens galaxy, 19°, is reasonably cons istent with the measured position angle of 26° ± 5°, 
which is typical for gravitational lenses (IKeeton et al.lll998l : iKoopmans et al.ll2006l ). Based 
on the Einst ein radius of the m odel, the Faber- Jackson relations for gravitational lenses 
measured by iRusin et al.l (120031 ) predict a lens galaxy apparent magnitude of R = 19.24 
(/ = 18.57) that agrees well with the measurement of i? = 19.32 ± 0.16 (J = 18.43 ± 0.25). 

Good fits t o quadruple lenses generally require both the ellipticity of the lens and an 
external shear (IKeeton et al.l 119971 ) . so for our second model we added a shear to the SIE 



model (the shear amplitude 7 and its position angle 9^). In this case, the best fit model 
overfits the data {Xrcd — 0.062 for dof = 1). Although the resul ting large amplitude of th e 
external shear is not inconsistent with N-body simulations (e.g.. [Holder fc Schechteiil2003l ). 
the lens galaxy position angle is misaligned with respect to the observations. Moreover, the 
large ellipticity and external shear cross almost perpendicularly, which indicates that the 
model may not be realistic. If we add weak constraints to match the axis ratio and position 
angle of the SIE component to the visible galaxy (e = 0.28 ± 0.15 and 6e = 26° ± 10°), 



then we obtain a good fit (x 



red 



1.2 for dof=3) with a small external shear (7 = 0.02). 
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This suggests that external shear is not important for fitting the image positions despite the 
possible existence of nearby groups. 

In all these models, the predicted flux ratios differ from the observed flux ratios (see 
Table [2] and Figure [6]) . This remains the case if we a dd the flux ratios an d their measurement 



errors as model constraints. The cusp relation (see lKeeton et al.ll2003l ) provides a means of 
determining whether the flux ratios of the three cusp images are consistent with any smooth 
mass distribution. In this case we find -Rcusp = 0.13 and (I/Re = 1.25 based on the i?-band 
flux ratios and the SIE model, which is marginally consistent with the range of distributions 
found for smooth lens models. The origin of the problem is presumably substructure in the 
gravitational potential of the lens due to either microlensing by the stars or sub-halos, since 
the observed flux ratios show no signiflcant wavelength dependence. We note, however, that 
evidence for anomalous flux ratios is not conclusive, mainly because of large astrometric and 
photometric errors. 



5. Summary 

We report the discovery of the quadruply imaged quasar SDSS J125H-2935. The lensing 
hypothesis is confirmed by the facts that i) the SDSS spectrum of the system shows the 
emission lines of a quasar at = 0.802 and the absorption lines of a galaxy at zi = 0.410, ii) 
the UH8k and Optic images confirm that the system consists of 4 blue point-like components 
and an extended object whose color is consistent with an elliptical galaxy at z ^ 0.4, iii) 
the geometry of the system is that of a typical cusp lens and can be well reproduced by a 
SIE model with reasonable parameter values, and iv) the luminosity of the lens galaxy is 
consistent with the expected luminosity from the Faber- Jac kson relation. Th is system is the 



second lowest redshift lensed quasar after RXJ1131-1231 (ISluse et al.ll2003l ) at Zs = 0.66. 
The fiux ratios of the three cusp images of the lens show a modest fiux ratio anomaly, whose 
origins could be better constrained with higher resolution images under better observing 
conditions. 
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Fig. 1. — SDSS i-band image of the SDSS J1251+2935 field (l'/2 seeing, 54 sec expo- 
sure). The pixel scale is 0'.'396 pixel"^, North is up and East is left. An expanded view of 
SDSS J1251-I-2935 is shown in the inset at lower left. Small numbers beside galaxies are 
R — I colors from the Optic images. 
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Fig. 2.— SDSS spectrum of SDSS J1251+2935 with a resolution of ~ 1800. The quasar 
emission hues (Mg II, [O II], and H7) redshifted to z = 0.802 are marked by the dashed hues 
and the galaxy absorption lines (Ca H&K, G-band, Mg, and Na) redshifted to z = 0.410 are 
marked by the dotted lines. The galaxy absorption lines indicate the existence of a bright 
lensing galaxy. 
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Fig. 3. — Optic /-band image of SDSS J1251-I-2935. The upper-left panel is the original data. 
The galaxy-like extended object in the upper-right panel is the residual after subtracting only 
4 point-like components, and vice versa in the lower- left (overplotted with contours). There 
are no residuals after subtracting a galaxy and 4 PSFs, as shown in the lower-right panel. 
The images and contours are scaled by the square-root of the counts. 
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Fig. 4.— Original Optic (BRI) and UH8k {V) images of SDSS J1251+2935 (upper panels) 
and the residuals after subtracting the best models (lower panels) . The best models consist 
of 4 PSFs for the Optic S-band and UH8k K-band, and 4 PSFs plus a galaxy for Optic R 
and /-bands. The pixel scales of Optic and UH8k are 0'.'274 pixel"^ (2x2 binned) and 0'.'232 
pixel"^, respectively. Although there are some residuals in the l/-band image due to the 
galaxy, its flux is too small to be measured accurately. 
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Fig. 5. — Color-color diagram {B — R and R~I) of the 4 point-l i ke com ponents. The arrow 
upper-left indicates the extinction direction for a lCardelli et al.l (119891 ) extinction law with 
Rv = 3.1 and A{B - V) = 0.3. 
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Fig. 6. — Flux ratios of the 4 point-like components. The flux ratios predicted by the 
SIE model (without flux ratio constraints) and the SIE plus external shear model (with the 
i?-band flux ratio constraints and the weak constraints on the ellipticity and the position 
angle of the lens galaxy) in §1] are plotted as dashed and dotted lines, respectively. Note that 
the latter model fits the observation with Xred = 2.5. 
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Fig. 7. — Predicted and observed positions of components A-D and G. The filled circles 
with error bars represents the observed positions, and the gray filled triangles represent the 
positions predicted by the SIE model. The critical curve {black line) and the caustics {gray 
line) are also plotted. The gray filled triangle on the caustics is the predicted source plane 
position of the quasar. 



Table 1. ASTROMETRY AND PHOTOMETRY OF SDSS J1251+2935 



Object 


x[arcsec]'^ 


?/[arcsec]^ 


5^ 






A 


-1.40 ±0.03 


-1.00 ±0.02 


20.67 ±0.22 


20.05 ±0.23 


19.41 ±0.10 


B 


-1.67 ±0.05 


-0.65 ±0.02 


21.06 ±0.31 


20.07 ±0.25 


19.72 ±0.27 


C 


-1.77 ±0.07 


0.04 ±0.04 


22.26 ±0.30 


21.48 ±0.28 


20.38 ±0.14 


D 


0.00 ±0.04 


0.00 ±0.04 


21.97 ±0.16 


21.33 ±0.16 


20.51 ±0.14 


G 


-1.04 ±0.09 


-0.23 ±0.07 




19.32 ±0.16 


18.43 ±0.25 



'^Measured in the Optic i?-band image using GALFIT. The positive directions of 
X and y are West and North, respectively. Errors indicate the dispersions from 6 
different PSF templates. 

''Measured in the Optic images using GALFIT. The errors are the dispersions 
from 6 different PSF templates, and they do not include the absolute calibration 
uncertainties. The magnitudes are calibrated using the standard star PG0918±029 
(iLandoltlEggj ). 



Table 2. SDSS J1251+2935: MASS MODELS 



Model 


Data^ 




2 c 
Aflux 


it!E[arcsec] 


e 


^e[deg] 


7 


0,[deg] 


comments 


SIE 


pos 


3.3/3 


(20) 


0.88 


0.19 


19 






bad flux 


SIE+shcar 


pos 


0.062/1 


(18) 


0.79 


0.66 


10 


0.21 


-80 


large misalignment 


SIE+shear 


pos+shape 


3.7/3 


(16) 


0.87 


0.21 


25 


0.018 


-8.7 


bad flux 


SIE 


pos+flux 


17/6 


8.8 


0.80 


0.46 


18 






poor fitting 


SIE+shear 


pos+flux 


8.9/4 


1.1 


0.76 


0.67 


53 


0.28 


-32 


large misalignment 


SIE+shear 


pos+flux+shape 


15/6 


6.3 


0.83 


0.38 


34 


0.074 


-15 


poor fitting 



Note. — Results of various mass models constrained by _R-band data. The position angles are measured East of 
North. The time delay between images A and D in the SIE model is AtAD ~ 17 /i~^day. 

^Data used to constrain the models; pos: positions of the 4 images and the galaxy, fiux: fiuxes of the 4 images, and 
shape: weak constraints on the ellipticity and position angle of the lens galaxy (e = 0.28 ± 0.15 and 9 — 26° ± 10°). 

''Total and the degree of freedom. 



■^Contribution of fluxes to the x^. The values in parentheses are not included in Xtot- 



